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Abstract. Observations are inconsistent with a homogeneous and iso- 
tropic universe with ordinary matter and gravity. The universe is far 
' from exact homogeneity and isotropy at late times, and the effect of 

the non-linear structures has to be quantified before concluding that 
^ ■ new physics is needed. We explain how structure formation can lead 

' to accelerated expansion, and discuss a semi-realistic model where the 

C/3 timescale of the change in the expansion rate emerges from the physics 

of structure formation. 



J> 1 Introduction 

^ 1.1 Three assumptions and a factor of 2 

m • 

CN| ■ The early universe, at least from one second onwards, is well described by a model 
I , which is homogeneous and isotropic (up to linear perturbations) and has only ordi- 
' nary matter and gravity. Here ordinary matter refers to matter with non-negative 
OO , pressure and ordinary gravity is based on the four-dimensional Einstcin-Hilbert 
' action. However, predictions of such models disagree with observations of the 
^ , late-time universe by a factor of two or so. The angular diameter distance to the 
' last scattering surface at a redshift of 1100 measured from the cosmic microwave 
, background (CMB) is about 1.7 times smaller than in the spatially flat homoge- 
■ neous and isotropic model dominated by pressureless matter. (Assuming that the 
primordial perturbations are adiabatic and that their spectrum is a power-law, 
and normalising to the same value of the Hubble parameter today.) Observations 
of type la supernovae and baryon acoustic oscillations at redshifts of order one 
and below are consistent with the CMB, and indicate that the discrepancy arises 
at late times, when the age of the universe is some billions of years. 
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In the context of the homogeneous and isotropic Friedmann- Robertson- Walker 
(FRW) models, there is a unique correspondence between the expansion rate and 
the distance scale (given the spatial curvature), and longer distances imply faster 
expansion. In terms of the expansion rate, compared to expectations, the Hubble 
parameter is a factor of 2 larger, given the observed energy density of matter 
(r^mo = 87rGN/3mo/(3i?o) ~ 0.25), or a factor of 1.2-1.5 larger, given the age of 
the universe (i?o^o ~ 0.8-1 instead of Hoto — 2/3). 

The factor 2 is observationally beyond reasonable doubt, so at least one of the 
three assumptions of the model is wrong. Either there is exotic matter with neg- 
ative pressure ( "dark energy" ) , or standard general relativity is not valid at large 
distances, or the universe is not homogeneous and isotropic, or some combination 
of the three. 

There is no evidence for exotic matter or modified gravity apart from ob- 
servations of the distance scale (and the related expansion rate). For example, 
modifications of gravity have not been observed in the solar system (apart from 
possibly the Pioneer anomaly) (Reynaud 2SQS^, and dark energy has not been 
measured in the laboratory. The situation is different from that of dark mat- 
ter, for which there are several independent lines of evidence, such as the CMB 
peak structure, gravitational lensing, the rotation curves of galaxies, the motions 
of clusters, the early formation of structures and so on. In contrast, the various 
measurements of the distance scale (or the expansion rate) are probing the same 
physical observable using different tracers. An effect which increases the expansion 
rate (and correspondingly the distance scale) at late times can account for all of 
the observations. 

Unlike the assumption of ordinary matter and gravity, the assumption that the 
universe is homogeneous and isotropy is known to be violated at late times, when 
non-linear structure formation has started. It is therefore necessary to quantify 
the effect of the breakdown of homogeneity and isotropy on the expansion rate 
before concluding that the observations call for new physics. It is suggestive that 
non- linear structures become significant at late times, when large departures from 
the predictions of the homogeneous and isotropic models with ordinary matter 
and gravity are observed. In contrast, models of exotic matter or modified gravity 
typically have no natural scale corresponding to billions of years, which is known 
as "the coincidence problem" . 

In section 2 we explain how inhomogeneity and/or anisotropy can affect the 
expansion rate. In section 3 we discuss a toy model which helps to understand 
the physics of acceleration due to structure formation. In section 4 we discuss a 
semi-realistic model of structures, and demonstrate how the correct timescale of 
some billions of years emerges from the physics of structure formation. In section 
5 we conclude with a summary. This is an overview of work reported in (Rasanen 
I2006al I2006bl I2008al l2008bD. 
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2 Backreaction 

2. 1 The fitting problem 

The effect of inhomogeneity and/or anisotropy on the average expansion rate is 
known as backreaction. To be specific, consider a manifold which one can split 
into space and time. Take a fixed spatial domain and consider average quantities 
such as the energy density and the expansion rate in this domain. If the domain is 
inhomogeneous and/or anisotropic, the time evolution of these average quantities 
will be different than in a homogeneous and isotropic space, even if the initial 
conditions are the same. 

One might expect the effects of inhomogeneities and/or anisotropics to be small 
in the real universe, since the universe looks homogeneous and isotropic when 
viewed on sufficiently large scales. However, one has to distinguish between ex- 
act and statistical homogeneity and isotropy. The early universe is locally close to 
exact homogeneity and isotropy, but after perturbations grow non-linear and struc- 
ture formation starts, the universe is only statistically homogeneous and isotropic. 
In other words, the local evolution is non-linear, but the distribution of the non- 
linear regions is homogeneous and isotropic, with a homogeneity scale which is 
today of the order 100 Mpc (Hogg I2004p . A space which is only statistically 
homogeneous and isotropic expands differently than a locally homogeneous and 
isotropic space, since averaging and time evolution do not commute. This non- 
commutativity was first discussed in detail by George Ellis (|1984D . and he termed 
the issue of finding the model which best fits the average of the locally complex 
behaviour "the fitting problem" . 

2.2 The local equations 

It is possible to give some rather general results on the average expansion rate. 
Let us consider a general spacetime where the matter can be treated as dust, that 
is, a pressureless ideal fiuid. After structure formation starts, this approximation 
does not hold (Pueblas & Scoccimarro 120081) . but the effects of velocity disper- 
sion related to small scale structures are not expected to be important for our 
discussion. The geometry is determined by the Einstein equation 



where Gaf3 is the Einstein tensor, Gn is Newton's constant, Tq,^ is the energy- 
momentum tensor, p is the energy density and u" is the velocity of observers 
comoving with the dust. 

We are interested in the average evolution of the geometry determined by ()2.1|) . 
In curved spacetime, only scalars can be straightforwardly averaged, so we will look 
at the scalar part of (|2.ip . By taking the trace we obtain one scalar equation, and 
projecting with and u^V^ gives two more (Buchert I1999P : 



(2.1) 




47rGNp- 2a2 + 2^2 



(2.2) 
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= 8^GNP-^(')i? + a2-c^2 (2.3) 

P + Op = 0, (2.4) 

where a dot stands for derivative with respect to proper time t measured by ob- 
servers comoving with the dust, dt = u^da, is the expansion rate of the local 
volume element, = \(j"'^(Tap > is the scalar built from the shear tensor CTq/j, 
= ^uj°'^ijJap > is the scalar built from the vorticity tensor ujafj, and ^^^R is 
the Ricci scalar on the tangent space orthogonal to the fluid flow. The equations 
are exact, and valid for arbitrary large variations in density, expansion rate and 
other physical quantities. The Hamiltonian constraint (|2.3p is the analogue of the 
Friedmann equation, and gives the local expansion rate, and (12. 4p expresses the 
conservation of mass. The Raychaudhuri equation (|2.2p gives the local accelera- 
tion. In the FRW case, we would have (j"^ = up — Q and = 3i/, where H is the 
Hubble parameter. 

We take the vorticity to be zero. While rotation is generated on small scales 
(Pueblas & Scoccimarro 2008), we expect that vorticity, like velocity dispersion, 
is not important for the physics we discuss. Given lo^ — 0, we see from (|2.2I) that 
the local acceleration is always negative (or at most zero). One might be tempted 
to conclude that since the local expansion decelerates everywhere at all times, the 
average expansion must also decelerate. However, this conclusion is false. 

2.3 The Buchert equations 

Let us see how the average expansion rate behaves by averaging the exact equations 
((^ - (P^ . This was first done by Thomas Buchert p999p . The average of a scalar 
quantity / is simply its integral over the spatial hypersurface of constant proper 
time, divided by the volume. 



/d3,V(3)5(t,x) ' ^'-'^ 

where '-^-'^ is the determinant of the metric on the hypersurface of constant proper 
time. 

Averaging (|2.2p - (|2.4p . we obtain the Buchert equations 

3- = _47rG'N(p) + Q (2.6) 
a 

3^ = 8^Gn(p) - - is (2.7) 

9,(p)+3-(p)=0, (2.8) 
a 

where the backreaction variable Q contains the effect of inhomogeneity and anisotropy: 
Q^\{{e^)-{9Y)^2{a^) , (2.9) 
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and we have defined a scale factor a{t) such that the volume of the spatial hyper- 
surface is proportional to a{t)^, 



where a has been normalised to unity at time ig: which we take to be today. 
As gives the expansion rate of the volume, this definition of a is equivalent to 
3a/a = (6). We will also use the notation H = a/a. 

The Buchert equations (|2.6p - (|2.8p differ from the FRW equations by the pres- 
ence of the backreaction variable Q and the related fact that the spatial curvature 
is not necessarily proportional to a~^. Because of this extra degree of freedom, the 
equations are not closed, but they provide a constraint on the expansion history. 
The backreaction variable Q in (|2.9p has two parts. The shear is also present in 
the local equations (|2.2p - (|2.3p . and acts to decelerate the expansion. The other 
term is the variance of the expansion rate, which has no counterpart in the lo- 
cal equations, and may in this sense be called emergent. This term expresses 
the non-commutativity of time evolution and averaging, and makes it possible 
for the average equations to have qualitatively different behaviour than the local 
equations. The contribution of the variance to the acceleration is always positive 
(except in the case of homogeneous expansion, when it is zero) . If the variance is 
large enough compared to the contribution of the shear and the energy density, 
the average expansion rate accelerates, as (j2.6p indicates, even though locally the 
expansion decelerates everywhere. 

3 Demonstrating acceleration 

3.1 A two-region toy model 

The reason why the average expansion rate can increase even though the local 
expansion decelerates is simple. In an inhomogeneous space, different regions ex- 
pand at different rates, and the volume of faster expanding regions increases more 
rapidly, by definition. Therefore the fraction of the volume which is expanding 
faster increases, and the average expansion rate can rise. 

Structure formation in the universe starts from a density distribution which is 
very smooth, with only small local variations. In a simplified picture, overdense 
regions slow down as their density contrast grows, and eventually turn around 
and collapse to form stable structures. Underdense regions in turn become ever 
emptier, and their expansion rate increases (relative to the mean). 

In order to clarify the physics of the acceleration due to inhomogeneity and/or 
anisotropy, we will consider a simple toy model of structure formation presented 
in (Rasanen I2006al I2006bp . before discussing a semi- realistic model in the next 
section. We consider two spherical regions, one overdense and one underdense, 
with scale factors denoted by oi and 02, respectively. In the Newtonian regime, 
a spherically symmetric region expands like a FRW universe with the same mean 
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density, regardless of the density profile. This is known as the spherical collapse 
model in the overdense case. 

We take the underdense region, which models a cosmological void, to be com- 
pletely empty, so it expands like Oi cx t. The evolution of the overdense re- 
gion, which models the formation of a structure such as a cluster, is given by 
a2 oc 1 — cos (j), t (X (f) — sin 0, where the parameter 4> is called the development 
angle. The value = corresponds to the big bang singularity, from which the 
overdense region expands until = tt, when it turns around and starts collapsing. 
The region shrinks to zero size at = 27r. In studies of structure formation, the 
collapse is usually taken to stabilise at (p = 3tt/2 due to vorticity and velocity 
dispersion, and we will also follow the evolution only up to that point. 

The total volume in the toy model is a'^ = af -|- a^. The average expansion rate 
and acceleration are 

H = ^^Hi + ^^H2=viHi+V2H2 (3.1) 

1 'd Hi Hi (Hi ~ H2f , , 

= -;^- = ]^^'i9i + ;^«2<Z2-2z;ii;2 , (3.2) 

where we have expressed the acceleration in terms of the deceleration parameter q 
(note that negative q corresponds to positive acceleration) . The average expansion 
rate is the volume- weighted average of the expansion rates Hi and H2, and it is 
therefore bounded from above by the fastest local expansion rate. In contrast, 
the acceleration is not simply the average of the local accelerations qi and q2. 
In addition to the first two terms in (|3.2p . which are the local values, there is 
a third term, which is non-positive, and contributes towards acceleration. This 
term depends on the difference of the expansion rates of the two regions, and 
corresponds to the variance part of the backreaction variable Q in (|2.6p . 

The toy model has one free parameter, the relative size of the two regions at 
some time. We fix this by setting the deceleration parameter q to be equal to 
that of the spatially flat ACDM FRW model when JIa = 0.7. In figure [T] (a) we 
plot 9 as a function of the development angle 4>. In addition to the toy model, we 
also show the ACDM model for comparison. This is done to help understand the 
toy model, and does not imply that the toy model should be taken seriously in a 
quantitative sense. 

The ACDM model starts matter-dominated, with q= 1/2. As vacuum energy 
becomes important, the model decelerates less and then crosses over to accel- 
eration. Asymptotically, q approaches —1 from above as the Hubble parameter 
approaches a constant. The toy model also starts with the FRW matter-dominated 
behaviour. As the overdense region slows down, the expansion decelerates more. 
Because the expansion of the overdense region is slower, the underdense region 
eventually takes over and comes to dominate the expansion rate. The expansion 
then decelerates less, and eventually accelerates; in fact the acceleration is stronger 
than in the ACDM model. 

The acceleration is not due to regions speeding up locally, but due to the slower 
region being less represented in the average. This is particularly easy to understand 
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(a) (b) 

Fig. 1. The evolution of the toy model as a function of the development angle <j). (a): 
The deceleration parameter q in the toy model (blue, solid) and in the ACDM model 
(red, dash-dot), (b): The Hubble parameter multiplied by time, Ht, in the toy model 
(blue, solid) and in the ACDM model (red, dash-dot). 



after the overdense region has started collapsing at (/) = tt. Then the contribution 
of the overdense region V2H2 to (|3.ip is negative, but the magnitude shrinks rapidly 
as V2 decreases, so it is transparent that the expansion rate increases. Note that 
while there is an upper bound on the expansion rate, there is no lower bound on 
the collapse rate. Therefore, the acceleration can be arbitrarily rapid, and q can 
even reach minus infinity in a finite time. (This simply means that the collapsing 
region is so dominant that vanishes in the denominator of q.) This is in contrast 
to FRW models, where q> ~1 unless the null energy condition (or the modified 
gravity equivalent) is violated. 

Figure [T] (b) shows the Hubble parameter multiplied by time as a function of 
the development angle <j). This is the same information as in figure [1] (a), but 
plotted in terms of the first derivative of the scale factor instead of the second 
derivative. In the ACDM model, Ht starts from 2/3 in the matter-dominated era 
and increases monotonically without bound as H approaches a constant. In the toy 
model, Ht falls as the overdense region slows down, then rises as the underdense 
region takes over, approaching Ht = 1 from below. The Hubble parameter in the 
toy model is at all times smaller than in the ACDM model, and because H is 
bounded from above by the fastest local expansion rate, Ht cannot exceed unity. 
This bound is valid also in realistic models: as long as the matter can be treated 
as dust and vorticity can be neglected, inhomogeneous and/or anisotropic models 
predict Ht < 1, in contrast to FRW models with exotic matter or modified gravity. 
Observationally, Ht k, 0.8-1.0 today (see section 5.1 of (Rasanen l2008ap ). 

Note that in an inhomogeneous and/or anisotropic dust model, whether the 
expansion accelerates simply depends on how rapidly faster expanding regions 
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catch up with slower ones: roughly speaking, how steeply the Ht curve rises. 
This is why the variance in Q (12. 9p contributes to acceleration in the average 
Raychaudhuri equation (|2.6p . With a larger variance, the difference between slow 
and fast regions is bigger, and faster regions take over more rapidly. 

4 Towards reality 

4.1 A semi-realistic model of bacl<reaction 

The toy model shows that acceleration due to structure formation is possible, and 
helps to understand what this means physically. This possibility has also been 
demonstrated with the exact spherically symmetric dust solution, the Lemaitre- 
Tolman-Bondi model (Chuang et al. 120051 Kai et al. 120061 Paranjape & Singh 



The issue is then whether this possibility is realised in the universe. The non- 
linear evolution of structures is too complex to follow exactly. However, because 
the universe is statistically homogeneous and isotropic, knowing some statistical 
properties is enough for evaluating the average expansion rate. In terms of the 
Buchert equations (I2.6p - (|2.8p . the average expansion rate is completely determined 
once the variance and the shear are known. The average expansion rate is also 
determined if we know which fraction of the universe is in which state of expansion 
or collapse at each time. We will now discuss a semi-realistic model which does this 
by extending the two fixed regions of the toy model to a continuous distribution 
of regions which evolves in time fRasanen I2008al l2008b[) . 

The starting point is the spatially flat matter-dominated FRW model, with a 
linear Gaussian field of density fluctuations. Structures are identified with spher- 
ical peaks (and troughs) of the density field smoothed on some scale, following 
Bardeen ([1986]). The peak number density as a function of the smoothing scale 
and peak height can be determined analytically. We keep the smoothing threshold 
fixed such that cr(t, i?) = 1, where a is the root mean square density contrast, t is 
time and R is the smoothing scale. Non- linear structures typically form at cr « 1, 
so R corresponds to the size of the typical largest structures, and grows in time. 
The smoothing is just a simplified treatment of the complex stabilisation and evo- 
lution of structures in the process of hierarchical structure formation. Since the 
peaks are spherical and isolated, and they are individually assumed to be in the 
Newtonian regime, their expansion rate is the same as that of a FRW universe 
with the same density, as in the toy model. The fraction of volume which is not 
in peaks is taken to expand like the spatially flat matter-dominated FRW model. 

The average expansion rate can be written as 



where vs (or, more accurately, dvg/dS) is the fraction of the volume in regions with 
linear density contrast 5 and expansion rate Hs{t). The correspondence between 
S and Hg is given by the spherical evolution model (i.e. the FRW evolution), and 
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Fig. 2. The BBKS (blue, solid) and BDG (red, dashed) transfer functions as a function 

of fc/fceq. 

the distribution of regions vs{t) is given by the peak statistics, which is determined 
by the power spectrum of the Gaussian density field. The power spectrum consists 
of the primordial spectrum, determined in the early universe by some process such 
as inflation, and the transfer function, which describes the evolution of the modes 
at later times. 

We take a scale-invariant primordial spectrum and the cold dark matter trans- 
fer function. For the latter, we consider two different approximations in order to 
show the uncertainty in the calculation. The BBKS transfer function (Bardeen 
I1986P is a fit to numerical calculations (we take a baryon fraction of 0.2), and the 
BDG form introduced by Bonvin et al. (|2006p is a simple analytically tractable 
function with the correct qualitative features. In figure [2] we show the suppression 
of the amplitude of modes as a function of k/keq, the wavenumber divided by the 
wavenumber of the modes which enter the horizon at matter-radiation equality. 
The amplitude of modes which enter the horizon during the radiation-dominated 
phase is suppressed, while modes entering during matter domination retain their 
original amplitude. The only scale present in the model is the one associated with 
the turnover near matter-radiation equality. With the transfer function fixed, there 
are no free parameters: the expansion history H(t) given by ()4.ip is completely 
fixed. 

In figure [3] we show Ht as a function of r = fccq^, the smoothing scale divided 
by the matter-radiation equality scale. As in the toy model, Ht « 2/3 at early 
times. As time goes on, ever larger non-linear structures form, and they take up 
a larger fraction of the volume. The expansion rate grows (relative to the FRW 
value) slowly, until there is a rapid rise and saturation, roughly at the scale of 
matter-radiation equality. It is clear that after r — I, when the perturbations 
which correspond to the matter-radiation equality scale collapse, Ht must settle 
to a constant, since the transfer function is essentially unity, and there is no scale 
in the system anymore. 

The matter-radiation equality scale is k^^ w IB.Tct;"^ Mpc « 100 Mpc, using 
the (somewhat model-dependent) value uj„i = 0.14 (Komatsu l2008p . Observation- 
ally, a{t, i?) ~ 1 today on scales somewhat smaller than 8 Mpc, so Rq w 10 
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(a) (b) 

Fig. 4. The expansion rate Ht as a function of time for (a) the BBKS transfer function 
and (b) the BDG transfer function. 



Mpc. Therefore the present day happens to be located around r = OT in the plots 
- exactly in the transition region. 

In figure [4] we show Ht as a function of time (in billions of years) instead of 
the smoothing scale. The timescale is set by the matter-radiation equality time 
teq and the amplitude A of the primordial density perturbations, t oc A^'^/^toq. 
The matter-radiation equality time is teq ~ lOOOw"^ years ^ 50 000 years for 
Wni — 0.14. The primordial amplitude of the density perturbations inferred from 
the observations depends on the damping of the perturbations after deviations 
from the matter-dominated FRW evolution become important. We take the same 
amplitude as in the ACDM model with I^a = 0.7 today, A = 4.6 x 10"^ (Dodclson 
I2003p . Ignoring damping by taking the matter-dominated FRW value A — 1.% x 
10~^ would make the timescale longer by a factor of 3.8. The slope of the Ht 
curve is less steep as a function of time than as a function of r, because the size 
of structures grows more slowly at late times. 

The timing of the rise of Ht is remarkably close to the observed acceleration 
era, and the change is of the right order of magnitude, 15-25%. Nothing related 
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to the present time has been used as an input in the model, the timing comes from 
the slope of the transfer function and the primordial amplitude. 

As noted when discussing the toy model, acceleration is a quantitative question 
related to the slope of the Ht curve. In the present case, while the expansion rate 
increases relative to the FRW case, the rise is not sufficiently rapid to correspond 
to acceleration. This is related to the fact that, unlike in the toy model, the 
overdense regions play almost no role, and the evolution of Ht can be understood 
entirely in terms of the underdense voids. At early times, voids take up only a 
small fraction of the volume, and Ht rises smoothly as they become prominent. 
In order for the change to be more drastic, the expansion rate should slow down 
due to overdense regions before the voids take over, as in the toy model. 

The model involves a number of approximations, and it should not be trusted to 
more than an order of magnitude. For example, we have not taken into account in 
the statistics that peaks located inside larger peaks should not be double counted, 
and we should also exclude troughs which are inside peaks, while the reverse is not 
true (overdense regions can be located inside voids) , which breaks the symmetry of 
the initial Gaussian density field (Sheth & van de Weveaert I2003P . Also, treating 
the peaks as isolated may not a good approximation once the peak density is high, 
collapsing structures are typically highly non-spherical, even the more accurate 
BBKS transfer function has an error of about 30% and so on. A more realistic 
calculation with quantified errors is necessary to estimate the effect of structures 
on the expansion rate in detail. 

5 Conclusions 

5.1 Summary 

Observations are inconsistent with a homogeneous and isotropic universe with 
ordinary matter and gravity. FRW models do not include the effect of non-linear 
structures on the expansion rate. The Buchcrt equations, where the effect of 
structures is taken into account, show that the average expansion of a space can 
accelerate due to inhomogeneity and/or anisotropy. 

With a simple toy model involving one overdense and one underdense region, 
we have demonstrated how acceleration is possible. First the expansion rate slows 
down due to the overdense region, and when the underdense region takes over, the 
change is sufficiently rapid that the expansion accelerates. The physical interpre- 
tation is simple: the fraction of volume in the faster expanding region grows, so 
the average expansion rate rises. 

The important question is whether this happens in the real universe. We have 
discussed a semi- realistic model which involves an evolving distribution of regions, 
with statistics given by the peak number density of the initial Gaussian density 
field. In the model, the expansion rate rises at some billions of years by about 15- 
25%, the right order of magnitude. The physical reason is that the volume of the 
universe becomes dominated by voids. The model has no acceleration, because 
overdense regions are never prominent. This is not crucial, since the model is 
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probably not correct beyond an order of magnitude, due to various approximations, 
such as the spherical symmetry of the structures and having equal amounts of 
mass in the overdense and underdense regions. It is remarkable that the correct 
timcscalc emerges from the physics of structure formation, essentially from the 
primordial amplitude of fluctuations and the time of matter-radiation equality 
imprinted on the cold dark matter transfer function. However, a more realistic 
calculation is needed to determine the details of the effect of structure formation 
on the expansion rate. 

One important aspect which we have not discussed here is the relationship 
between the average expansion rate and observations of distance and redshift. 
Physically, it is clear that as the volume expands more rapidly, distances will grow 
longer. However, the quantitative relationship between the expansion rate and the 
distance scale remains to be determined (this is work in progress). 

That structure formation could account for the observations of the late-time 
distance scale and expansion rate is a plausible possibility. Until the effect is quan- 
tified in detail, we do not know whether new physics is needed or the observations 
can be understood in terms of a complex realisation of general relativity. 
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